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Abstract 

Diffusion of symmetric diblock copolymer chains in macroscopically oriented 
lamellar block copolymers are studied in a molecular dynamics simulation. 
Results for diffusion constant both parallel JDm and perpendicular D± to the 
lamellar planes are compared to results in the disordered one phase region. 
For diblocks of length N = 40 (l.lN e ) and 100 (3N e ), where N e is the en- 
tanglement length of a homopolymer melt at the same density D\\ is nearly 
independent of \, while D± is strongly suppressed as \ is increased. These re- 
sults are in agreement with theoretical predictions based on the Rouse model. 
The isotropic diffusion constant of quenched disordered systems is approxi- 
mated fairly well by (2-Dii + D±_)/3 of the corresponding lamellar system. 
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The question whether lamellar ordering modifies self diffusion in a symmetric diblock 
melt has been the subject of several recent theoretical and experimental studies. Barrat and 
Fredrickson |]J argued that for chains obeying Rouse dynamics, the diffusion within a lamel- 
lar plane is unaffected by the lamellar structure, while perpendicular diffusion is retarded by 
an amount determined by the degree of segregation. Measurements on quenched, unoriented 
samples [[H§ indicated that there is no discontinuity for the diffusion constant at the order- 
disorder transition (ODT). The samples used in these studies were polycrystalline, without 
macroscopic ordering, so that the separate components of the diffusion coefficient could not 
be determined. Such a separation was achieved [|J for entangled chains by macroscopically 
orienting the sample using oscillatory shear. Measurements on such mono domain samples 
showed less anisotropy than predicted by the Rouse model. For well entangled chains, both 
parallel and perpendicular diffusion have to involve translation of the A segments of the 
chains into the B domains, so that both components are reduced by the lamellar structure, 
though not necessarily to the same extent. Such an activated reptation mechanism in which 
the diffusion constant D decays as exp(— xN) was observed by Lodge and Dalvi [[5] for 
highly entangled melts of poly(ethylenepropylene)-poly(ethylethylene) diblock copolymers. 
The difference between the diffusion constants parallel D\\ and perpendicular D±_ to the 
interface was found to decrease as the ODT is approached. Block-retraction mechanism, 
similar to arm-retraction known from star polymers, is also found to play a role in the 
diblock diffusion and also gives rise to an exponential decay of D with increasing N. 

Computer simulations of diblock copolymer systems have been mostly restricted to study- 
ing the statics and thermodynamics of such systems ||. The only simulation which studied 
diffusion in a diblock system is by Pan and Shaffer 0, who studied a symmetric diblock in 
the strong segration regime (xN = 45) for several chain lengths N using the bond fluctu- 
ation method on a lattice. They found that D\\ was roughly two times smaller than D for 
a homopolymer melt and the critical chain length at which entanglements became impor- 
tant was comparable in the two systems. No perpendicular diffusion was observed for this 
value of xN. Since these simulations were done on a lattice, there was no a priori way to 
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determine the equilibrium lamellar spacing d\. 

Here we study by continuum space molecular dynamics simulations the diffusion con- 
stants of perfectly oriented lamellar systems both in the plane of the lamella and in the 
perpendicular direction, for chains of length 40 and 100, corresponding to about 1.1 and 
3N e respectively. N e ~ 35 is the entanglement length for a homopolymer melt j5|, at the 
same average density (p ~ 0.85cr~ 3 ). d\ is determined from a constant pressure simulation. 
It is then fixed and D» and D± are determined at constant volume. The results are compared 
to simulations in the isotropic disordered phase. 

We use a coarse grained bead-spring model similar to that used in our earlier studies 
of polymer melts and networks || and tethered chains ||. Each chain consists of N beads 
(monomers), connected to form a linear chain. For diblocks, the model is generalized to 
two types of polymer species, say A and B, which are connected in a block of fN beads 
of type A connected to a block of (1 — f)N beads of type B. Here we consider only 
symmetric diblocks, / = 1/2. The interaction potential Ujj(r) between two beads of types 
I, J = {A, B} separated by a distance r is taken as the repulsive part of a Lennard Jones 
6 : 12 potential, 
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where the cutoff distance r c = 2&aij. Here, ejj and 07 j are, respectively, parameters fixing 
the energy and length scale for monomers of type / and J. Adjacent monomers along the 
chains are coupled through an anharmonic FENE potential, chosen to eliminate unphysical 
bond crossings or breaking |8]||. All the interactions are short ranged and the model is very 
efficient computationally. For symmetric copolymers, we set gaa — £bb — £ab/c — 1 —I - c, 
and oaa — &bb — g ab — a - The latter choice simply means that all monomers have the same 
interaction range and thus the two species have the same molar volume. All of our results 
are reported in terms of a and e, with r = a{mje) x l 2 for the time scale, where m is the mass 
of a monomer. The energy parameters are a special case of (e = 1) e = €ab — \ {.£aa + £bb), 



used in simple lattice models such as in Flory-Huggins theory. In ref. ||10|| ) we showed that 
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simply introducing this difference in the repulsive interaction strength between like and 
unlike monomers was sufficient to drive phase separation. By varying e, while leaving the 
temperature T constant, one can induce a phase transition from a homogeneous blend to 
a system with two segregated coexisting phases [e c = 3A0(5)kBT/N] [[HJ. This approach 
also has the advantage of overcoming a problem inherent in experimental studies: namely, 
the dependence of the monomeric friction constant on the temperature. It is somewhat 
reminiscent of recent elegant experiments by Gehlsen et al. in which the phase separation 



of binary mixtures is studied as a function of the difference in deuterium content between 
two otherwise identical polymer species. 

We study M chains of length N by solving Newton's equations of motion using a velocity- 
Verlet algorithm [12] with a time step At = 0.013r [§,[L(J. Two different starting states were 



used: disordered and lamellar. In the former case, we studied cubic systems at density 
p = 0.85a 3 [f|. The lamellar state was constructed as described in ref. [pTD|| . Since di 



depends on e and is not known a priori, we first ran a constant pressure simulation JED] at 
a P = 5ecr~ 3 in which the dimensions of the simulation cell, parallel L» and perpendicular 
to the lamellar plane varied independently |]T3|| for each value of e. di adjusts itself 



rather rapidly (a few hundred thousand time steps). With increasing e, L± increases and 
the cross sectional area L| decreases, such that the overall density remains fairly constant. 
The equilibrium lamellar spacing increases rapidly with e and then saturates as shown in 
Fig. [I]. The results scale with N as predicted by self-consistent field theory [f[4-|16l . For large 



e the chains are highly stretched and the number of AB contacts is small. Thus, further 
increase in the repulsive energy between the two species does not strongly affect di. In this 
regime x is no longer proportional to e as often incorrectly assumed in many previous lattice 
simulations of diblock copolymer systems |||7|. 

After di is determined, we fixed the volume to study the dynamics in the lamellar phase. 
The results presented here are for M = 1600 chains of length N = 10, M — 800 chains of 
N = 20 and 40 and M = 200 chains of iV = 100. The iV = 40 and 100 systems were studied 
in both the disordered and lamellar regions of phase space, while N = 10 and 20 systems 



were only simulated in the disordered phase. The lamellar systems contained 100 chains 
per layer for N = 40 and 50 chains per layer for iV = 100, giving L± = 4di for N = 40 
and 2di for N = 100. From our study of the order parameter in the ordered phase as a 
function of e, we estimate that e O Dr(40) = 0.85 ± 0.05 and e O Dr(100) = 0.28 ± 0.03 |10 . 
To study the diffusion in both phases, long runs were necessary to obtain good statistics. 
Typical runs for N = 40 were of length (5.2 — 10.0) x 10 4 r while for iV = 100, they were of 
order (3.9 — 6.5) x 10 5 r. In addition, lamellar systems for N = 200 and N = 400 were also 
simulated for visualization of the chain motion. However, these simulations were too large 
to run long enough to obtain the diffusive constants. These results will be presented in a 
later paper. 

We follow the motion of the chains' center of mass through gs(t), whose components are 
defined by 



with a = x, y, or z. Here r cma (t) are the components of the center of mass of a chain and 
(• • •) denotes both an ensemble average over all the chains in the system and an average 
over different starting states. At sufficiently long times, g^ a {t) exhibits diffusive motion with 
gz a — > 2D a t. In the disordered phase, the three diffusion constants are identical and equal 
to D, the isotropic diffuion constant. In the lamellar phase, the motion in the xy plane 
(plane of the lamella) is expected to be faster than the motion between the lamellar planes 
(z direction). The parallel diffusion constant Dm = (D x + D y )/2, and D^ = D z . 

The behavior of {g^ x {t) + g^ y (t))/2 for some of the N = 40 systems is shown in Fig. |2|. 
For the isotropic systems (homopolymers and the disordered diblocks), the quantity shown 
is gz(t)/3. For times longer than about 10 4 r, all the systems show diffusive behavior. The 
same quantity for a homopolymer melt and a lamellar system of chain length N = 100 is 
shown in Fig. [| The transition to the diffusive regime in the lamellar system takes place 
around t « 10 4 r, somewhat later than in the homopolymer melt. This delay in the onset of 
the diffusive motion is common to all lamellar systems studied. The diffusion constants were 
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calculated from the slope of (g 3x (t)+g 3y (t))/2 in the diffusive regime. Perpendicular diffusion 
is observed only for systems sufficiently close to the ODT. For large e, <?3 Z (t) saturates at 
about the square of the thickness of the interface region, within which the junction points 
of the chains are localized. 

Chains of length 40 (l.liV e ) are clearly in the Rouse regime and are expected to obey 
Rouse dynamics. Barrat and Fredrickson modelled the self diffusion of a Rouse chain in 
a lamellar system as tracer difusion of such a chain in a periodic potential varying in one 
dimension, perpendicular to the lamella planes 0. According to this study, diffusion in 
the plane is not affected by the periodic potential, while diffusion perpendicular to the 
lamella is slowed down by an amount dependent upon the magnitude of the concentration 
fluctuations. Near the ODT, D± is expected to decrease by less than 40% compared to D\\, 
while in the strongly segregated phases it should decay exponentially with xiV []J . To check 
these predictions, we plot in Fig. |] both D\\ and D± as a function of e for lamellar systems 
as well as the isotropic D for disordered systems. D» in the lamellar N = 40 systems is 
only very weakly dependent on e. At the ODT, D\\t/<j 2 ~ 9.5 x 10~ 4 , while in the strong 
segregration limit (e = 9.0, not shown in Fig. [|), D\\r/a 2 decreases only to 8.0 x 10~ 4 . For a 
homopolymer melt (e = 0), the isotropic diffusion constant is about 1.5 x 10~ 3 o~ 2 /r, so that 
the diffusion constant in the isotropic phase close to ODT is lower by about 50% compared 
to the diffusion constant in the homopolymer melt. This reduction results from the increased 
friction from the local motion of the monomers. D±, on the other hand, decays rapidly with 
e. While near the ODT the ratio D±/D\\ ~ 0.8, it is reduced to about 0.03 for e = 1.4. The 
reduction of about 20% near ODT is consistent with the Barrat-Fredrickson prediction. For 
higher values of e, we could not observe any perpendicular diffusion during the time scale of 
our simulations. 

The four non-zero values of D± for N = 40 can be fitted nicely with an exponential 
decay of the form D±/D a = 0.57exp[— 0.16iV(e — Iodt)], where D Q is the value the diffusion 
constant in the homopolymer melt. To make a quantitative comparison with the theory of 
Barrat and Fredrickson we note that they predict that when x4>o ~ Xodt, where ipo is 
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the amplitude of the density variation in the lamella, D± is reduced by a factor of about 10. 
For our system at e = 1.0 > £odt, = 0.83, so that eip Q = 0.83 ~ Iodt- For this system, 
D± is about 1/3 of the ODT value or about 3 times less suppression than predicted. 

The situation for the N = 100 chains is qualitatively similar. This length corresponds to 
about 3N e , and entanglements are not expected to play a major role. In a homopolymer of 
N = 100 chains, crossover to reptation motion is barely observed [fy. For stretched chains, 
N e might even increase. The systems near the ODT all have essentially the same D\\ (about 
2.1 x 10 - V 2 /t). The reduction in D at the ODT compared to the homopolymer melt is 
about 40%, comparable to that seen for N = 40. For very high values of e, the reduction in 
D\\ is comparable to that for N = 40 chains. In the most strongly stretched system studied 
(e = 9.0, not shown in Fig. f|), D\\ is about 0.75 of Dii near ODT, while the same ratio is 
about 0.85 for iV = 40. The perpendicular motion is significantly reduced in comparison to 
the parallel one. For e = 0.3 (close to the ODT), we found D± = 7 x 10 _5 <7 2 /t, or about 
1/3 of the corresponding Dn. We note that the suppression of D± near ODT is higher than 
for the N = 40 diblocks and the Barrat-Fredrickson prediction. However, for these longer 
chains, a more precise determination of the ODT is needed before a reliable comparison can 
be made. For larger values of e, no measurable perpendicular diffusion was found. For the 
system with e = 0.5, only two chains succeeded in translating their junction points from one 
interface to another during the whole course of the simulation, with no measurable D±. For 
larger e, no chain exchange occurred between the interfaces. 

It is also interesting to note that for quenched disordered systems (those with e > £odt), 
the isotropic D ~ (2-Dii + D±)/3 of the corresponding lamellar system jnj. Fredrickson 
18] had suggested that for small anisotropy, the diffusion constant of a quenched system 



should be slightly larger than (2D\\ + D±)/3. We do not see an evidence of the predicted 
enhancement. However, this prediction is only for systems very close to the ODT. Since the 
exact ODT of our model system is not known to a high accuracy, we cannot rule out this 
result. 

We also studied the diffusion in the disordered regime for shorter, unentangled chains 



of length N = 10 and 20. For iV = 10 and 20, the isotropic diffusion constant scales 
D ~ N^ 1 f{e/eoDT)i where f(x) is a decreasing funciton of x. A comparision of these 
results with the theoretical predicitions of Leibig and Fredrickson [[19] for tracer diffusion in 
the disordered regime will be presented in a later paper. 

Our results demonstrate that for lamellar systems of short chains (of length N = 40 
and 100), self diffusion in the lamellar plane near the ODT is essentially unaffected by the 
amount of chain stretching. This is in accordance with the Rouse-like motion of chains. 
Perpendicular diffusion is strongly suppressed by the lamellar structure. Our results are 
qualitatively consistent with the Barrat-Fredrickon theory, although the amount of suppres- 
sion is somewhat less than predicted. In the limit of strong stretching, some reduction in D\\ 
is observed, possibly indicating the onset of entanglement effects. Further studies on long 
chains are needed to study the effect of entanglements. 

We thank T. P. Lodge for helpful discussions. This work was supported in part by the 
BMBF project "Computer Simulations of Complex Materials". 
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FIG. 2. Mean squared diplacement of the center of mass parallel to the lamellar plane, 
[(d3x(t) + 93y(t))/2] versus t/r for TV = 40 systems: lamellar with e = 0.9 (A) and 3.0 (o), 
homopolymer melt (■), isotropic near the ODT, e = 0.8(D), and at e = 3.0(«) . 
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FIG. 3. (gsx(t) + 93y(t))/2 versus t/r for a homopolymer melt (o) and a lamellar system with 
e = 1.0 (□) for N = 100. The lines have a slope of 1, characteristic of the diffusive regime. 
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FIG. 4. Diffusion coefficients as a function of the excess AB interaction parameter, e: D\\ (o) 
and D ± (□) for iV = 40 lamella, isotropic D (A) for disordered TV = 40 melt, D\\ (o) for iV = 100 
lamella and isotropic D (y) for disordered N = 100 melt. The lines for the N = 40 systems are 
only a guide to the eye. 
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